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Abstract 


The Cuernavaca aquifer is the main source of permanent fresh water that 
supports economic development and supply to the population of the 
capital of the state of Morelos. During the last thirty years in the south of 
the Cuernavaca aquifer nitrates have increased degrading groundwater 
quality. The identification of the origin of the nitrate is essential to 
implement public policies that minimize the contamination. Taking into 
account the characteristics of the region such as soil type and geology, 
population growth, economic activities such as agriculture, and 
recreational activities including golf, four possible sources of nitrate are 


proposed in the study area: 1) Natural sources (soil organic nitrogen), 2) 
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infiltration of wastewater, 3) infiltration of agricultural water, and 4) 
infiltration of water used in the irrigation of golf courses. In the southern 
portion of the Cuernavaca aquifer, the geochemistry and environmental 
isotopes of water (9180, 9%H, 33H) indicate that the main nitrate 
contaminant source is sewage disposal. The data also indicate that the 
arsenic detected in the wells comes from the deepest strata of the aquifer. 
The results also show correlation between  nitrate and  uranium 
concentration. The data obtained allowed to determine the meteoric water 


line of the Cuernavaca Aquifer. 


Keywords: Cuernavaca aquifer, nitrate, sewage water, geochemistry, 
environmental isotopes of water, arsenic, uranium, meteoric water line of 


Cuernavaca aquifer. 


Resumen 


El acuífero Cuernavaca es la principal fuente de agua dulce permanente 
que sustenta el desarrollo económico y el suministro a la población de la 
capital del estado de Morelos. En los últimos 30 años, los nitratos en el 
sur del acuífero Cuernavaca se han incrementado, degradando la calidad 
del agua subterránea. La identificación del origen del nitrato es 
fundamental para implementar políticas públicas que controlen y 
reviertan dicha contaminación. Con base en las características propias de 
la región, como son tipo de suelo y geología, incremento acelerado de la 
población, y actividades industriales, agrícolas y recreativas, entre ellas 
la práctica del golf, se plantean cuatro posibles fuentes de nitrato en la 
zona de estudio: 1) nitrato del suelo, 2) infiltración de aguas residuales, 


3) fertilizantes provenientes de zonas agrícolas, y 4) fertilizantes 
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provenientes de campos de golf. En la porción sur del acuífero 
Cuernavaca, los datos geoquímicos y los isótopos ambientales del agua 
(9180, 32H, 9H) indican que la contaminación por nitrato proviene 
principalmente de aguas residuales. Asimismo, indican que el flujo 
profundo aporta arsénico a algunos aprovechamientos subterráneos. Los 
resultados también muestran correlación entre las concentraciones del 
nitrato y del uranio en el acuífero. El registro periódico de los datos 
isotópicos permitió también determinar la línea de agua meteórica de la 


zona del acuífero Cuernavaca. 


Palabras clave: acuífero Cuernavaca, nitrato, aguas residuales, 
geoquímica, isótopos ambientales del agua, arsénico, uranio, línea de 


agua meteórica del acuífero Cuernavaca. 


Received: 30/12/2021 
Accepted: 29/11/2022 
Published Online: 13/12/2022 


Introduction 


Nitrate (NO3”) is one of the main causes of groundwater contamination 
(Zhao, 2015; Canter, 1997). Various sources of nitrate contribute to this 
contamination, including agricultural and livestock activities, and 
wastewater (Haller, McCarthy, O'Brien, Riehle, 8 Stuhldreher, 2013; 


Viers, Liptzin, Rosenstock, Jensen, 8 Hollander, 2012). 
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The consumption of drinking water with high concentrations of 
nitrate is a risk factor for several types of cancer, including gastric, 
colorectal and brain tumors (CDPH, 2013). The World Health Organization 
(WHO) established a maximum permissible limit of 50 mg/l for nitrates 
(WHO, 2007) to prevent damage to the health of the population. For 
Mexico, the permissible limit (Modification to the Official Mexican 
Standard NOM-127-SSA1-1994) (SSA, 1994) for drinking water is 10 
mg/l of nitrogen (as N-NO3”), which is equivalent to 44.3 mg/l of nitrate. 


Aquifers contaminated with nitrates have been reported in Mexico, 
both from agricultural activities and wastewater (Pacheco 8 Cabrera, 
1997; Cardona, 2004; Torres-Martínez et a/., 2021). In the study carried 
out by Suarez-Barragán and Ureno-Luna (1988) about the quality of 
groundwater in the Cuernavaca Valley, the presence of nitrates was 
detected in wells in some municipalities; among them, the Tetecalita well 
in Emiliano Zapata with 6.2 mg/l, a spring in Temixco with 4.43 mg/l, a 
spring in Cuernavaca with 6.64 mg/l, the Melchor Ocampo well with 5.75 
mg/l and the Hotel Casino de la Selva well with 5.85 mg/l, both in 
downtown Cuernavaca. By 1995 (Conagua, 2002), nitrate concentrations 
in the western portion of the Cuernavaca aquifer (municipalities of 


Temixco and Acatlipa), already reached values of 20 mg/l (Figure 1). 
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About 89 Mm of wastewater is generated annually in Morelos, of 
which 66 % is discharged without treatment to septic tanks, absorption 
wells, or directly into streams, ravines and rivers (Conagua, 2010). 
Discharges are mainly public-urban with 46 % and services with 21 %. 
The municipal wastewater flow is approximately 4 585 l/s, of which 2 820 


l/s are collected and 959 l/s are treated (Conagua, 2010). 


58 % of Morelos' surface area is used for agricultural and pasture 
production. The entity has 204,441 hectares dedicated to the agricultural 
sector (INEGI, 1991). In the study area, the productive crop areas are 
located in Jiutepec, in the ejidos of Acapantzingo, Chapultepec and 
Atlacomulco; to the northeast, in the limits of the municipalities of 
Cuernavaca and Tepoztlán; to the south, in the ejidos of Temixco, 
Acatlipa, Zapata and Tezoyuca (POZCI, 2009). The main crops grown in this 
region are tomatoes, green tomatoes, corn, rice and ornamental horticulture 
(Sagarpa, 2003). Pesticides, herbicides and fertilizers used in agricultural 
areas, when they infiltrate into the subsoil, can contaminate groundwater 
(Conagua, 2010). 


Four golf courses are located in the study area. Some of the 
nutrients applied to the grass are (Del-Campo-Delgado, 2016): a) 
Nitrogen, an essential element for grass growth and giving green color; 
b) phosphorus, related to root development, and c) potassium, which acts 
at many levels of plant metabolism. When water and nutrients are applied 


in excess, fertilizers infiltrate in the form of leachates. 


Nitrates in the southern Cuernavaca aquifer have increased over the 
last 30 years, degrading groundwater quality. The identification of the 


origin of nitrate is fundamental to implement public policies to control and 
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revert such contamination. 


The objective of this study was to determine nitrate concentrations 
in the southern portion of the Cuernavaca aquifer, and identify the 
sources that produce nitrate. To do so, hydrogeological tools, 
groundwater chemistry, and environmental water isotopes (9180, 3*H, 


93H) were used. 


Study area 


The Cuernavaca aquifer is located in the physiographic province of the 
Mexican Volcanic Belt, subprovince Balsas-Mexcala Basin (Raisz, 1964), 
between parallels 18% 40' and 19% 10' North latitude, and 99% 00' and 99% 
20" West longitude (Figure 2). 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


o MI) Check for updates 
OPEN (Y ACCESS 


Tecnología y 


CienciaszAgua 


59150" 


E 


Figure 2. Location of Cuernavaca aquifer. 


The geomorphological limits of the Cuernavaca aquifer are: to the 
north the Zempoala, Ajusco and Chichinautzin mountains ranges, to the 
south the Xochicalco mountain range, to the west the extension of the 


Zempoala mountain range and to the east the Tepoztlán-Tlaltizapán 


AE 024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


mountain range (Conagua, 2020). The Zempoala (3 650 masl) and 
Chichinautzin (3 690 masl) mountain ranges are the highest elevations 
(Conagua, 2020). 


According to mean annual temperatures, the Cuernavaca aquifer 
zone has four types of climates (Conagua, 2002): Semi-cold (5-12 *C), 
temperate (12-18 *C), semi-warm (18-22 *C) and warm (> 22 *C). The 
semi-cold and temperate climates are located in the northern portion, 
while the semi-warm and warm climates are located in the central and 


southern portions, respectively. 


The orography and relief of the terrain (Conagua 2002) influence 
rainfall that affects the Cuernavaca aquifer. In the upper portion of the 
region, precipitation is in the order of 1 500 to 1 100 mm per year, while 
in the southern portion it is in the order of 950 to 850 mm per year. In 


general, precipitation is more abundant from June to October. 


In the northern portion, various types of temperate and cold climate 
forests predominate, including: Pine forest, oyamel forest, cedar-tascate 
forest, oak forest, pine-oak forest and mountain mesophyll forest 
(Secretaría de Desarrollo Sustentable de Morelos, 2009). In the central 
and southern portions, low deciduous forest, low subtropical scrub and 
pasture predominate (Secretaría de Desarrollo Sustentable de Morelos, 
2009). Agricultural crops are also a type of non-natural vegetation 
present in the plain zone. 


The stratigraphic column of the Cuernavaca aquifer is composed of 
thirteen types of rock units, which range from Lower Cretaceous to recent 
(Conagua, 2020). These units are: Xochicalco Formation, Morelos 


Formation, Cuautla Formation, Mexcala Formation, Balsas Group, 
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Tilzapotla Rhyolite, Colotepec Granite, Tepoztlan Formation, Zempoala 
Andesite, Cuernavaca Formation, Chichinautzin Formation, Continental 


Clastic Deposits and Quaternary Alluvial Deposits. 


Hydrogeology 


Based on geomorphological and geohydrological aspects, it is possible to 
identify three hydrogeological units in the Cuernavaca aquifer (Figure 3) 
(Conagua, 2002; Morales-Casique, Guinzberg-Belmont, 8  Ortega- 
Guerrero, 2016; Ortega, Garcia, 8: Pelayo, 2003). 


Sierra 
Chichinautzin 


Spring Sierra 


Spring Palo Bolero Xochicalco 
San Ramón 
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The first unit corresponds to a free aquifer hosted in fractured 
basaltic igneous rocks of the Chichinautzin formation, and continental 
clastic deposits of the Cuernavaca formation. The second unit, defined as 
an aquitard, is associated with low permeability materials of the Balsas 
Group and the Mexcala Formation. It functions as the basement of the 
free aquifer, and is located at average depths of 100 to 200m in the 
middle and southern portion of the aquifer. The third unit is a semi- 
confined regional aquifer, hosted in the limestones of the Morelos and 
Cuautla formations, and is the basis of the geohydrologic sequence of the 


Cuernavaca aquifer. 


The preferential recharge zones are located in the upper parts of 
the study area. Sierra Zempoala in the northwest, made up of fractured 
volcanic rocks from the Tertiary, and Sierra Chichinautzin in the north, 
made up of fractured and cavernous basalts from the Quaternary (Jaimes- 
Palomera et al., 1989). 


Anthropogenic sources of pollution 


The main streams and ravines that drain wastewater from the city of 
Cuernavaca, in a north-south direction, are El Pollo ravine in the western 
portion, and Barranca Analco in the central and eastern portions of the 
municipality of Cuernavaca, which flow north of the municipal capital of 
Temixco into the Apatlaco River. The northeastern part of Cuernavaca 
discharges into the Puente Blanco ravine, which collects wastewater from 
Jiutepec and Emiliano Zapata, where it is called Arroyo Dulce, and finally 


flows into the Apatlaco River south of the town of Xochitepec. In the 
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easternmost portion of the valley drains the Las Fuentes ravine, into 
which the Huauchiles and Las Fuentes springs discharge, as well as 
wastewater from Tejalpa, Progreso, Las Fuentes and El Pochotal, finally 
reaching the agricultural areas of southern Jiutepec and Emiliano Zapata 
(Figure 4) (IMTA, 2012). 
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Methodology 


Rainfall 


Rainwater samples were collected from 2016 to 2022, using RS2 PALMEX 
rainwater collectors, which allow collecting accumulated samples without 
evaporation or isotopic fractionation (IAEA, 2002). Samples were filtered 
using cellulose membranes with 0.45 um pore size, collected in high- 
density polyethylene bottles with conical and polyseal insert caps, and 


preserved at 4 *C. 


The collectors are located in the following towns: Tezoyuca, 
Temixco, Cuernavaca, Huitzilac, Tres Marías and Tres Cumbres. These 
collectors are part of the National Network of Isotopes in Precipitation 
(Known as RENIP), which covers the 15 physiographic provinces that 
make up the national territory. The RENIP is operated by the department 
of Hydrology of the Mexican Institute of Water Technology (IMTA). 


Stable isotope analyses were conducted at the isotope hydrology 
laboratory of IMTA using laser water isotope analyzers Picarro L2110-i 
and Los Gatos Research LWIA-45EP. The analytical precision of the 
measurements was + 0.5%o (Vs. VSMOW) for 92H and + 0.13%o (Vs. 
VSMOW) for 9180. 
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Groundwater 


Water samples were collected from 16 groundwater developments during 


the months of January and February 2014 (Figure 5). 
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Samples were taken at the discharge train of the wells. Prior to 
measuring field parameters and collecting samples, the entire volume of 
water in the well annulus was renewed by pumped extraction, in order to 


sample and measure representative water from the aquifer. 


Samples were filtered using cellulose membranes with 0.45 um pore 
size and collected in high-density polyethylene bottles with conical and 


polyseal insert caps. 


Samples for anions and stable isotopes were preserved at 4 *C. 
Samples for metals and cations were acidified to pH < 2 with ultrahigh 


purity nitric acid. 


In order to reduce as much as possible the modification of the 
prevailing physicochemical properties of the water in the aquifer, isolation 
cells were used for the measurement of field parameters. These devices 
are connected with hoses to valves located in the discharge train of the 
well. This eliminates the contact of the water with the atmosphere, 
avoiding modifications due to the release of gases occluded in the water, 


such as carbon dioxide and/or hydrogen sulfide. 


Field parameters 


The parameters measured in the field were: hydrogen potential (pH), 
electrical conductivity (EC), dissolved oxygen (DO), redox potential (Eh), 
water temperature and ambient temperature (T). The methodologies 


used are: 
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1. The measurement of hydrogen ion activity in groundwater was 
carried out by means of a potentiometer. The calibration of the 
equipment was performed using the three-point calibration 
technique, from buffer solutions of known pH. In this technique, the 
pH to be recorded must be located between these extreme values. 
Solutions of 4.0, 7.0 and 10.0 pH units were used in the calibration 
of the equipment. 

2. Electrical conductivity was measured by means of a portable 
conductivity meter with temperature compensation. Calibration of 
the equipment was performed using a standard potassium chloride 
(KCI) solution of 1 412 uS/cm at 25 *C, 

3. Measurement of dissolved oxygen in groundwater was carried out 
using a temperature-compensated, coated galvanic membrane 
sensor electrode. The dissolved oxygen electrode was calibrated 
according to the atmospheric conditions at the site, ¡.e. in equilibrium 
with the air at site elevation. 

4. The redox potential of groundwater was measured with a combined 
type electrode, consisting of a platinum filament and a silver-silver 
chloride reference electrode with 3 mol/l potassium  chloride 
reference electrolyte. The sensitive element where the electron 
transfer takes place is the platinum filament, so it was kept clean and 
polished to avoid erroneous measurements. A standard Zobell's 
solution, Eh of 428 mV (a mixture of potassium ferrocyanide and 
ferricyanide) was used to verify the operation of the electrode and 
adjust the field measurement. 

5. Ambient and water temperature were measured with electrodes 


previously calibrated in the laboratory and verified in the field. 
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Parameters measured in the laboratory 


The measurement of dissolved chemical constituents in the groundwater 
samples was carried out by ¡on chromatography for anions and by two 
instrumental techniques for cations, depending on their concentrations: 
plasma-induction coupled with optical emission spectroscopy (ICP/OES) 
and plasma-induction coupled with mass spectrometry (ICP/MS). The 
ICP/MS technique is approximately 1,000 times more sensitive than the 
ICP/OES technique. The analyses were performed by Activation 
Laboratories Ltd of Canada (ACTLABS). For analytical quality control of NO3 
(as N), the laboratory used a certified reference standard of 2.87 mg/l. In 
the control tests applied to the measurements, the analytical error was less 
than 5 %. 


Results 


Rainfall 


The arithmetic average and precipitation-weighted average values of 
the 010 (%o) and 32H (%o) are shown in Table 1. Precipitation- 
averaged values (Sp) were calculated as dp=2P¡*8;¡/Pr, where P¡ is the 
amount of rainfall, 3; is the isotopic composition of the sample (32H or 


09180), and PT is the amount of total precipitation. 
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Table 1. Values of the arithmetic average and precipitation-weighted 


average values of the 9180 (%o) and 32H (%o). 


a Ain Longitude | Altitude 9180 9H 9:80» 9*Hp 

O (masl) (Y%o0) | (Y%o) | (%o) (%o) 

Tres Cumbres | 19.061267 | -99.218317 3290 -10.84 | -70.46 | -12.18 | -81.22 
Huitzilac 19.032850  -99.268800 2600 -9,43 |-60.99 -10.71 | -71.38 
Tres Marías 19.051483  -99,246283 2813 -9,31 |-59.69 | -10.4 | -68.91 
Cuernavaca 18.942833 | -99,215139 1638 -8.6 |-57.83| -9.99 | -68.12 
Temixco 18.857528 -99.223917 1264 -7,45 -49.85 | -8.95 -61.11 
Tezoyuca 18.800383 — -99.196733 1194 -7,43 |-50.41 | -8.74 -60.05 


Figure 6 presents the precipitation-weighted average of the 3180, 


The most enriched values in heavy isotopes correspond to the site with 


lower altitude and warm temperature, while the most depleted values 


correspond to the area with higher altitude and cold temperature. 
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Figure 6. Average values of the 0180 in the collectors. 


With the isotopic data of the samples, plots of the ratio 9180 %o (Vs. 
VSMOW) versus 9H %o (Vs. VSMOW) were elaborated. Figure 7 shows 
the data for the collectors located in the Chichinautzin mountain range, 
and Figure 8 shows the data for the valley. Using a least-squares linear 
regression of the isotope data, the Meteoric Water Lines were determined 
which are: 0H = 8.15 010 + 17.86 (Tres Cumbres); 0%H = 8.15 0180 + 
16.24 (Tres Marías); 0%H = 8 9180 + 14,58 (Huitzilac); 92H = 7.7 0180 + 
8.14 (Cuernavaca); 9%H = 7.53 9180 + 6.28 (Temixco); 0%H = 7.42 9180 


+ 4.7 (Tezoyuca). 
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Figure 7. Meteoric water lines of the collectors in the Chichinautzin 


mountain range (Tres Cumbres, Tres Marías and Huitzilac). 


92024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


Tecnología y 


CienciaszAgua 


20 


Q Cuernavaca 


0%H =7.7*5180 + 8.14 
R? = 0.9844 


-20 O Temixco 
z 92H = 7.53*5180 + 6.28 
S 40 R? = 0,9944 
YN 
> 
S -60 
o 
-80 
-100 
-120 


a 1) Check for updates 
OPEN ACCESS 


LMM 


(O Tezoyuca 


9%H = 7.42*5180 + 4.7 
R?= 0.9921 


6 -4 -2 0 


5180 (%o ) VSMOW 


Figure 8. Meteoric water lines of the collectors in the valley 


(Cuernavaca, Temixco and Tezoyuca). 


Taking into account the data from all collectors (Figure 9), the 


equation of the Meteoric Water Line of the Cuernavaca Aquifer zone ¡s: 
9%H = 7.8 910 + 8.82, R?= 0.9841, n = 118. 
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Figure 9. Meteoric Water Line of the Cuernavaca Aquifer. 


When evaporation takes place in environments with relative 
humidity less than 100 %, the fractionation of +0 exceeds that of ?H, 
producing vapor with high deuterium content (excess deuterium) relative 
to 180 (Martinelli, Victoria, Sternberg, Ribeiro, 8 Moreira, 1996; Gat, 
2000). This effect, defined as d= d2H-8*0918%0, produces a deuterium 
excess (d) of 10%o at the Global Meteoric Water Line (GMWL) at relative 
humidity of 85 % (Merlivat 8, Jouzel, 1979; Clark 8 Fritz, 1997). 


The Meteoric Water Lines of the mountain range collectors (>2600 
masl) have intercept values (+17.86 %o, +16.24 %o, +14.58 %o) greater 
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than the GMWL value (+10 %o). The deuterium excess (d) values of these 


collectors are also greater than those of the GMWL (Figure 10). 
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Figure 10. Deuterium excess values (d) of rainfall from collectors 


located in the Chichinautzin mountain range. 


Vapor generated by evapotranspiration from forested areas has 
deuterium excess values greater than 10 (Galewsky et a/., 2016). When 
this vapor condenses, it produces rain with deuterium excess values 
greater than 10. Therefore, the high values of deuterium excess and 
rainfall intercept of Tres Cumbres, Tres Marías, and Huitzilac indicate that 


evapotranspiration is an important source of moisture in the Chichinautzin 
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mountain range. Gastmans et al. (2017) found a similar effect in Brazil. 
The Meteoric Water Lines of the Belo Horizonte, Rio de Janeiro, Campo 
Grande, Carolina, Betania and Rio Claro stations presented intercept 
values greater than +10 %o. They attributed this effect to the 
recirculation of vapor generated in the Amazon rainforest. Esquivel- 
Hernández et a/. (2019) also reports this effect in a tropical ecosystem in 


Ecuador. 


In arid atmospheres, when rainfall is light or just beginning, 
raindrops on their descent from the clouds to the ground undergo 
evaporation (Peng, Mayer, Harris, € Roy, 2007). This is known as 
secondary evaporation or evaporation under the cloud. Fractionation 
produces enrichment of the 3180 and 3H, and decreases the deuterium 
excess (Froehlich, Kralik, Rank, 8 Scheifinger, 2008; Stewart, 1975; 
Dansgaard, 1964). When precipitation is intense or long duration, the 


environment becomes saturated and secondary evaporation stops. 


In the case of the Cuernavaca valley, with two collectors located in 
very warm climate (Temixco and Tezoyuca), the lowest deuterium excess 
values correspond with the smallest precipitations, indicating secondary 
evaporation of light rain (Figure 11). With respect to the most abundant 
rainfall, their deuterium excess values are similar to the GMWL value (+10 
%o), indicating that these samples did not evaporate (Figure 11). 
Therefore, the Meteoritic Lines of the Cuernavaca valley collectors, have 
intercept values (+8.14 %o, +6.28 %o, +4.7 %o) lower than the GMWL 


value (+10 %o), due to secondary evaporation. 
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Figure 11. Values of deuterium excess (d) versus precipitation, of the 


collectors located in the valley of the Cuernavaca Aquifer. 


The altitudinal isotopic gradient is represented by the equation 
(Figure 12): 8180 = -0.0014 (Z)-7.1933 con R? = 0.8792 (Z in meters). 
Hence, the vertical gradient is -0.14 %o/100m. In Figure 12, the 3180 
values correspond to precipitation-weighted average. This vertical 
gradient is similar to that found by Scholl Gingerich and Tribble (2002) (- 
0.12 %o0/100m) for rainfall from Maui, Hawaii; and Otte et al. (2017) (- 
0.11 %o0/100m) for rainfall from Mt. Kilimanjaro, Tanzania. The altitude 


effect explains 87 % of the variation in isotopic composition (R? = 0.87). 
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Figure 12. Altitudinal isotopic gradient of the 0180, 


Groundwater 


Tables 2a, 2b, 2c and 2d shows the values of the physical-chemical 
parameters measured in the field and laboratory, of the samples from the 


south of the Cuernavaca aquifer. 


Tecnología y ciencias del agua, ISSN 2007-2422, 
Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


Table 2a. Physical-chemical parameters of the samples from the south 


of the Cuernavaca aquifer. 


Location Temperature 
ID Date td of the sample pH Ae E 
Xx Y Altitude (0) (ec) (4S/cm) | (mg/I) 

Well 1 30/01/2014 | 18%47.813 | 99%12.198 1186 21.2 19.4 6.85 641 0.30 

Well 2 30/01/2014 | 18%48.792 | 99%11.058 1217 28.5 20.0 6.86 268 5.60 

Well 3 31/01/2014 | 18%48.829 | 99%11.266 1220 23.3 22.7 7.07 620 7.30 

Well 4 31/01/2014 | 18%48.500 | 99%12.030 1202 23.8 19.4 7.45 506 0.60 

Well 5 04/02/2014 | 18%49.838 | 99%11.727 1219 23.7 24.0 6.75 902 6.50 

Well 6 04/02/2014 | 18%49.650 | 99%10.865 1241 24.1 22.8 6.88 674 4.30 

Well 7 05/02/2014 | 18%50.296 | 99%10.222 1241 25.5 22.7 6.81 831 3.90 

Well 8 05/02/2014 | 18%48.588 | 99*11.713 1203 31.1 19.8 7.40 281 0.80 

Well 9 06/02/2014 | 18%46.368 | 99%11.089 1138 23.9 20.6 6.80 279 5.20 

Well 10 06/02/2014 | 18%48.376 | 99%11.538 1280 22.9 21.1 7.13 697 1.50 

Well 11 06/02/2014 | 18%47.580 | 99%11.667 1182 26.7 19.5 7.29 777 0.90 
Spring San Ramón 08/02/2014 | 18%47,519 | 99*11.983 1189 27.0 20.5 6.80 1207 0.60 
Well 12 08/02/2014 | 18%46.893 | 99%12.576 1151 29.7 22.6 6.94 1270 4.30 
Spring Chihuahuita 18/02/2014 | 18%45.385 | 99%11.338 1096 30.7 20.4 6.69 392 4.80 
Well 13 19/02/2014 | 18%49.216 | 99%12.782 1232 27.6 22.4 6.97 542 3.50 
Spring Palo Bolero 18/09/2014 | 18%46.015 | 099014,312 1482 31.2 21.7 6.60 2240 1.08 
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Table 2b. Physical-chemical parameters of the samples from the south 


of the Cuernavaca aquifer (continued). 


1b E NO3 PO4 SO. Na Li Mg Al Si K Ca Ti v Mn 
ate 
(mg/l) | (mg/l) | (mg/l) | (mg/l) | (mg/l) | (mg/!) | (mg/!) | (mg/l) | (mg/l) | (mg/l) | (19/D | (u9g/l) | (mg/l) 
Well 1 30/01/2014 5.79 0.1530 145 13 0.028 27 0.020 25 2.2 100 <2 <2 0.0048 
Well 2 30/01/2014 16.71 0.0306 26 17 0.005 14 0.010 26 3.0 30 <i 14.8 0.0022 
Well 3 31/01/2014 17.11 0.1530 82 23 0.032 22 0.010 29 4.4 86 <i 11.7 0.0020 
Well 4 31/01/2014 9.95 0.1530 77 14 0.010 19 0.020 25 2.3 64 <2 4.2 0.0034 
Well 5 04/02/2014 42.08 0.1530 97 54 0.010 25 0.020 31 6.8 109 <2 4.6 0.0050 
Well 6 04/02/2014 28.02 0.1530 69 24 0.012 33 0.010 35 5.3 72 2.9 13.3 0.0005 
Well 7 05/02/2014 26.87 0.1530 102 29 0.013 38 0.010 34 5.3 96 2.8 12.6 0.0005 
Well 8 05/02/2014 11.54 0.1530 67 15 0.016 19 0.449 28 2% 57 24.1 12.6 0.0275 
Well 9 06/02/2014 16.84 0.5508 30 16 0.010 14 0.043 27 2.8 33 <2 10.8 0.0034 
Well 10 06/02/2014 14.19 3.0600 52 18 0.010 24 0.020 28 3.5 79 <2 3.6 0.0030 
Well 11 06/02/2014 6.19 0.1530 240 13 0.023 32 0.001 25 2.4 127 1.8 5.5 0.0055 
Spring San 
Ñ 08/02/2014 0.11 0.1530 429 16 0.054 49 0.020 24 2.8 230 <2 <2 0.0054 
Ramón 
Well 12 08/02/2014 9.50 0.1530 451 20 0.059 56 0.010 29 37 225 2.4 1.6 0.0005 
Spring 
18/02/2014 19.14 0.0306 33 14 0.005 13 0.112 23 2.5 32 2.6 15.7 0.0005 
Chihuahuita 
Well 13 19/02/2014 27.32 0.1530 45 22 0.031 36 0.010 37 2.6 34 2.7 14.1 0.0005 
Spring Palo 
sl 18/09/2014 | 0.0001 0.0001 1100 20 0.060 86 0.020 24 3.2 437 2.8 <2 0.0010 
olero 
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Table 2c. Physical-chemical parameters of the samples from the south 


of the Cuernavaca aquifer (continued). 


Error % 
Fe Co Cu As Rb Sr Mo Cs Ba TI Pb U Tritium 
ID Date ionic 
(mg/D) | (u9/D) | (mg/l) | (mg/D | (19/D) | (19/0 | (19/10) | (19/D | (19/D) | (u9/) | (mg/D) | (19/1) (UT) e 
Well 1 30/01/2014 0.1 < 0.1 0.002 0.05 5.7 584 10.0 1.34 35.6 0.800 0.00010 0.86 0.83 2.14 
Well 2 30/01/2014 0.05 0.070 0.002 0.0002 7.6 188 1.0 0.24 33.4 0.020 0.00200 0.76 0.79 2.87 
Well 3 31/01/2014 0.05 < 0.05 0.001 0.01 4.0 979 1.0 8.28 102.0 0.020 0.00005 1.91 2.17 10.67 
Well 4 31/01/2014 0.1 < 0.1 0.002 0.04 5.5 502 6.0 0.60 85.0 0.280 0.00010 1.25 0.76 6.86 
Well 5 04/02/2014 0.1 0.540 0.006 0.01 12.1 321 2.0 0.48 90.2 0.140 0.00200 2.31 2.16 2.32 
Well 6 04/02/2014 0.05 0.070 0.001 0.00 8.8 515 1.0 0.08 161.0 0.040 0.00005 2.76 1.35 2.23 
Well 7 05/02/2014 0.05 0.140 0.001 0.00 8.3 736 1.0 0.04 264.0 0.040 0.00005 4.76 2.11 0.03 
Well 8 05/02/2014 0.05 < 0.05 0.002 0.05 6.9 614 5.0 0.47 83.6 0.050 0.00200 1.43 0.12 0.08 
Well 9 06/02/2014 0.1 < 0.1 0.002 0.0003 7.1 224 <2 0.20 53.4 0.060 0.01000 0.95 0.75 4.67 
Well 10 06/02/2014 0.1 < 0.1 0.002 0.03 7.6 950 2.0 0.58 136.0 0.100 0.00010 2.16 1.53 1.03 
Well 11 06/02/2014 0.005 0.023 0.001 0.07 6.9 1430 10.0 1.25 29.2 0.440 0.00240 0.56 0.86 3.62 
Spring San 
Romá 08/02/2014 0.1 < 0.1 0.002 0.03 8.1 1550 2.0 2.46 21.2 0.160 0.00010 0.24 0.12 0.97 
amón 
Well 12 08/02/2014 0.05 0.050 0.001 0.03 11.8 1960 2.0 2.02 30.9 0.080 0.00005 0.25 0.86 1.11 
Spring 
18/02/2014 0.05 0.070 0.001 0.00 6.8 239 1.0 0.21 46.7 <0.01 | 0.00005 0.91 0.12 8.16 
Chihuahuita 
Well 13 19/02/2014 0.05 0.060 0.003 0.00 5.5 333 <1 0.14 39.0 <o0.01 | 0.00200 3.05 1.47 0.77 
Spring Palo 
Bj 18/09/2014 0.46 < 0.1 0.002 0.00 10.3 3180 <2 3.04 14.4 <0.02 | 0.00040 0.02 0.55 
olero 
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Table 2d. Physical-chemical parameters of the samples from the south 


of the Cuernavaca aquifer (end). 


Alkalinity Hardness 
Redox HCO3 SDT 510 5%H Deuterium F Cl NO» Br 
ID Date Total Total 
(mv) (mg/l!) | (mg/l) (%o) | (Y%o) | excess (d) | (mg/l) | (mg/l) | (mg/l) | (mg/l) 
(mg/l) (mg/l) 
Well 1 30/01/2014 152 212 259 462 361 -11.0 | -77.4 10.9 0.59 8 0.082 0.05 
Well 2 30/01/2014 477 121 148 193 131 -10.7 | -74.4 11.2 0.23 9 0.016 0.13 
Well 3 31/01/2014 349 168 205 446 305 -10.3 | -73.1 9.6 0.41 27 0.082 1.12 
Well 4 31/01/2014 453 211 257 364 240 -10.8 | -76.7 10.1 0.56 16 0.082 0.45 
Well 5 04/02/2014 354 256 312 649 375 -10.0 | -71.7 7.9 0.60 86 0.082 3.71 
Well 6 04/02/2014 370 232 283 485 315 -9.8 -70.6 8.2 0.23 39 0.082 1.34 
Well 7 05/02/2014 429 277 338 598 396 -10.1 | -71.9 8.8 0.21 58 0.082 2.28 
Well 8 05/02/2014 434 160 195 202 222 -10.8 | -75.7 10.7 0.88 20 0.082 0.79 
Well 9 06/02/2014 278 144 176 201 138 -10.6 | -74.4 10.4 0.26 12 0.016 0.17 
Well 10 06/02/2014 402 240 293 502 294 -10.3 | -73.8 8.3 0.05 36 0.033 1.21 
Well 11 06/02/2014 255 184 224 559 447 -11.0 | -77.8 10.6 0.74 8 0.082 0.05 
Spring San 
, 08/02/2014 36 328 400 869 775 -11.0 | -77.0 11.3 0.95 16 0.082 0.05 
Ramón 
Well 12 08/02/2014 320 323 394 914 791 -10.9 | -76.6 10.4 0.62 21 0.082 0.05 
Spring 
18/02/2014 425 144 176 282 133 -10.6 | -74.2 10.7 0.31 13 0.016 0.19 
Chihuahuita 
Well 13 19/02/2014 389 208 254 390 233 -10.4 | -73.4 9.4 0.28 17 0.082 0.05 
Spring Palo 
Bal 18/09/2014 120 307 375 1613 1448 -11.0 | -76.7 11.5 0.76 19 0.0001 150 
olero 


mg/!. 


permissible limit of 44.3 mg/l. The main anion is HCOz", with the following 


order of abundance HCO3"» SO4?>» CL:>NO3' (Figure 13). 
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Water temperature ranges from 19.4 to 24 *C; dissolved oxygen 
from 0.3 to 7.3 mg/l; pH from 6.6 to 7.45, and TDS from 193 to 1 613 


Nitrate concentrations range from O to 42 mg/l, below the 
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Figure 13. Anion composition of groundwater (% in meq/!). 


The main cation is Ca?*, with the following order of abundance Ca?*+ 
» Mg?* » Na* » K* (Figure 14). 
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Figure 14. Cation composition of groundwater (% en megq/!). 


Flow through the Chichinautzin formation produces Ca-Mg-HCO3 
type waters (Figure 15) (Chihuahuita spring, well 9 and well 2). Flow 
through the Cuernavaca and Balsas formations produces Ca-Mg-HCO3- 
S04 type waters (well 1, well 4 and well 8). Flow through the Morelos and 
Cuautla formations produces Ca-Mg-SO4-HCO3 type waters (Palo Bolero 
spring, San Ramon spring, well 11 and well 12). A group of wells (3, 5, 6, 
7, 10 and 13) have higher concentrations of chloride than the rest of the 


wells. 


2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


a 1) Check for updates 
OPEN ACCESS 


Tecnología y 

CienciaszAgua 
0 -99915' eo 
o o 
pu % 


Na+* Ci 
y E 
Mo sos 


1845" 
1845" 


Water types 


Calcium-bicarbonate 
Calcium-sulfate 
Calcium-magnesium 


ES Aquifer boundary 
Ó wells and springs 


-99915' 


Figure 15. Stiff diagrams of groundwater samples. 
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The concentration of nitrate (NO3") for 2014 showed levels lower 
than 5 mg/l in the eastern portion of the aquifer, where rocks of the 
Chichinautzin Formation predominate, composed mainly of fractured 
basalts (Figure 16). In contrast, the western portion showed the highest 
nitrate levels, with concentrations greater than 30 mg/l, where the 
Temixco-Acatlipa urban zone is located (Figure 16). Most of the urban 
zone of Cuernavaca, Jiutepec, and Emiliano Zapata, showed 


concentrations ranging from 10 to 20 mg/l (Figure 16). 
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Figure 16. Nitrate concentrations in the years 1995 and 2014. 
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It is important to note the higher concentrations of uranium in wells 

with high concentrations of nitrate (Figure 17). Ingestion of drinking water 
with high uranium concentrations can cause damage to the liver (Brugge 
8: Buchner, 2011) and bones due to alpha radiation (Kurttio et a/., 2005). 
The World Health Organization (WHO) established a maximum 
permissible limit of 15 ug/l for uranium (WHO, 2008). 
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Figure 17. Uranium versus nitrate. 


The oxidized states of uranium U(VI) are soluble in groundwater, 
while the reduced U(IV) species are immobile (Singh et a/., 2014; Nolan 
8: Weber, 2015; Paradis et al., 2016). Reduced U(IV) uranium species are 


very susceptible to being oxidized, when groundwater contains oxidizing 
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species such as nitrate (NO3”, y NO2”). Nolan and Weber (2015) found in 
two major aquifers in the United States, High Plains and Central Valley, a 
significant correlation between nitrate and uranium. They conclude that 
nitrate, the primary contaminant in these aquifers, should be considered 


as a factor inducing secondary uranium contamination. 


Stable isotopes 


The stable isotope data of oxygen and hydrogen in a scatter plot 9180 vs. 
92H allow us to identify and differentiate the types of water of each well 
(Figure 18). The graph identifies four groups of waters characterized by 
the altitude at which the effective recharge of groundwater takes place. 
The first group corresponds to wells that capture water of depleted 
isotopic composition, with 9180 values between -11.05%o and -11.02%o, 
and 0%H between -77.75%o and -76.71%o (wells 1 and 11, and San 
Ramón and Palo Bolero springs), due to the fact that their recharge zone 
is located at higher elevations than those corresponding to the rest of the 
wells. 
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Figure 18. Isotopic composition (318%0, 92H) of water from wells. 


The second group is characterized by 91%0 values between -10.69 
and -10.60 %o, and 02H between -75.6 and -74.1%o, which corresponds 
to the recharge zones with the lower elevation (wells 2, 9 and Chihuahuita 
spring). 

A third group, formed by wells 4, 8 and 12, which capture water 


with 0180 and 9H values intermediate between group one and two. 


The fourth group corresponds to waters that evolved isotopically 
from the previous groups, because they contain water components 


enriched in heavy isotopes due to intense evaporation. The isotopic values 
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are located around a line with a slope lower than the Cuernavaca Aquifer 
line (LMAC). This group is formed by wells 3, 5, 6, 7, 10 and 13. In the 
area where these wells are located there are extensive areas of irrigated 
crops, a golf course and several wastewater canals. The i¡sotopic 
compositions vary for 9180 between -10.35 and -9.84%o, and for 92H 
between -73.,4, and -70.6%o. This fourth group was classified as wells 


with evaporated water recharge. 


Figure 19 shows for the fourth group a proportional increase in 9180 
and Cl, evidencing recharge with water subjected to intense evaporation 


before infiltrating into the soil. 
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Figure 19. Values of 910 versus Cl. 
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Evaporation decreases the deuterium excess, while recycled 
moisture increases it (Froehlich et a/., 2008). Evaporation of a surface 
water body produces deuterium excess values lower than 10 (Victoria, 
Martinelli, Mortatti, 8 Richey, 1991). Therefore, Figure 20 indicates that 


water samples with higher NOz3 content contain evaporated water. 
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Figure 20. Values of NOz3 versus deuterium excess (d). 
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Tritium 
group of wells with evaporated water contain between 1.5 and 2.5 


tritium units (TU), indicating very short residence times, ¡.e., they receive 


modern water recharge (Figure 21). On the other hand, wells with non- 


evaporated water have tritium contents of less than 1 TU, indicating 


longer residence times. In the case of the San Ramón spring, the absence 


of tritium indicates water with residence times of more than thirty years 


after infiltrating the aquifer. 
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Figure 21. Values of 9180 versus tritium. 
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The water samples with the highest NO3 values contain between 1.5 


and 2.5 tritium units. This indicates very short subsurface residence times 
(Figure 22). 
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Figure 22. Values of NOz3 versus tritium. 
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Ionic relations 


The dissolution of gypsum produces a Ca?*/ SOa?- equal to 1: 


CaSO,:2H,0 + Ca?*+ SO,” + 2H,0 


Samples from the Palo Bolero spring, San Ramón spring, and wells 
11 and 12 are in the 1:1 equiline, indicating that gypsum dissolution is an 
important process in the chemical composition of these samples (Figure 
23). 
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Figure 23. Values of Ca versus SO. 
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The geochemical processes that control the chemical composition of 
water can be elucidated using the Gibbs diagram (Gibbs, 1970). In Figure 
24, total dissolved solids (TDS) are plotted as a function of the cation ratio 
(Na*/(Na* + Ca?**). Four groups can be distinguished in Figure 24. In the 
flow through the Morelos and Cuautla formations (red circles), the 
dominant process is the weathering of geologic material. In the flow 
through the Chichinautzin formation (green squares), the dominant 
process is rainfall. In the flow of the Cuernavaca formation (blue 
diamonds), its i¡onic composition is determined by weathering of the 
geologic material, but also by the effect of pluvial precipitation. In the 


fourth group (yellow triangles), the predominant process is evaporation. 
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Figure 24, Gibbs diagram: Cation ratio. 
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Linear regression between nitrate and chloride content yielded a 
value of R? = 0.531 for the evaporated samples (Figure 25). Correlation 
values greater than r?=0.35 indicate that the groundwater contains 
domestic or municipal wastewater (Ritter € Chirnside, 1984; Pacheco, 
Marín, Cabrera, Steinich, 8: Escolero, 2001). 


0.8 
0.7 


0.6 


NO; (meq/!) 
o o 
h u 


o 
w 


0.2 


O0.i 


0 0.5 1 1.5 2 2.5 3 
Cl (megq/!, 


Figure 25. Values of Cl versus NOz. 


The graph (molar ratio) of chloride versus nitrate normalized by 
chloride (NO3/C!), allows to elucidate the origin of the nitrate (Figure 26). 
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Samples fall in the blue box if the origin is agricultural, in the purple box 
if the nitrate is from manure, etc. The graph indicates that the nitrate in 
the evaporated samples comes from wastewater (green box) (Liu, Li, 
Lang, 8 Xiao, 2006; Guo, Yan, Wang, Xu, 8 Yang, 2020; Torres-Martínez 
et al., 2021). 
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Figure 26. Chloride (molar concentration) versus NO3/Cl (molar ratio). 
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The chloride and sodium values of the evaporated water samples 
(Figure 27) are similar to the values of septic tank effluent in the United 
States (Panno et al., 2002). 
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Figure 27. Chloride versus sodium. 


Samples that exceed the permissible limit for arsenic (NOM 127- 
SSA1-127 modified in 2000) (SSA, 1994) are those with low dissolved 


oxygen concentrations (Figure 28) and tritium values lower than 1 TU 
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(Figure 29), indicating that wells with high arsenic values are recharged 


by deep flows. 
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Figure 28. Values of Arsenic versus dissolved oxygen. 
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Figure 29. Arsenic versus tritio. 


Torres-Martínez et al. (2021) determined in the Comarca Lagunera, 
an agricultural region in northern Mexico, the sources of nitrate to 
groundwater: Manure (47 %), sewage (43 %), fertilizers (5 %), soil 
organic nitrogen (3 %), and atmospheric deposition (2 %). Their results 
indicate that even though the Comarca Lagunera is a very important 
agricultural area, synthetic fertilizers are not the main source of nitrate. 
Ren, Zhang, Wang and Wang (2021) conducted a study to determine the 
sources of nitrate in surface and groundwater in the Ye River basin of 


northern China. It is a predominantly agricultural area, whose populations 
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do not have sanitary sewers. For their study, they used stable isotopes of 
nitrogen and water, as well as a Bayesian mixing model. Their results 
indicate that the main contribution comes from wastewater and manure. 
They suggest building infrastructure for wastewater treatment. Pastén- 
Zapata, Ledesma-Ruiz, Harter, Ramírez and Mahlknecht (2013) 
conducted a study in a citrus-growing area in northeastern Mexico to 
determine the sources of nitrate in shallow groundwater. Based on 
geochemical evidence, they determined that the sources of nitrate in rural 
areas are the use of manure as fertilizer and wastewater, while in urban 


areas it is wastewater. 


Our results indicate that the nitrate present in the south of the 
Cuernavaca aquifer comes mainly from wastewater, which is consistent 


with sites with similar problems, such as those mentioned above. 


Nitrate in surface water and groundwater originates from multiple 
sources, and it is difficult to accurately distinguish these sources using 
hydrogeochemical methods alone. Nitrate isotopes have been shown to 
be a very useful tool for identifying nitrate sources (Torres-Martínez et 
al., 2021; Ren et a/l., 2021). In a new study being conducted in the 
Cuernavaca aquifer, hydrogeochemical methods as well as stable isotopes 


of water and nitrate are being used. 
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Conclusions 


This study presents the use of environmental water isotopes and 
geochemistry to determine the concentrations and origin of nitrate that 


degrade water quality in the southern Cuernavaca aquifer. 


The results show that the zones of the aquifer recharged with 
intensely evaporated water have the highest nitrate concentrations. The 
results also indicate that the nitrate comes from domestic or municipal 
wastewater, probably generated by untreated sewage discharges. None 
of the samples exceed the permissible limit of 44.3 mg/l of nitrate. 
However, the evolution of the last thirty years shows a significant increase 
in the area, which is why it is important to implement public policies that 
promote the construction and rehabilitation of infrastructure to collect 
wastewater and send it to existing treatment plants (La Gachupina and 
Emiliano Zapata WWTPs). 


The data obtained also revealed a correlation between nitrate and 
uranium. Even though all samples have uranium concentrations below the 
permissible limit, it is recommended to continue monitoring this 
parameter, since an increase in nitrate concentrations could be 
accompanied by an increase in uranium concentrations. Because uranium 
decay increases the radioactivity of groundwater, it is recommended to 
measure radioactivity to establish the current values and monitor its 


evolution over time. 


With respect to arsenic, the groundwater samples that exceed the 
permissible limit for arsenic have low dissolved oxygen and tritium 


concentrations, indicating that they come from the deep strata of the 
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aquifer. If groundwater extraction in the area increases in the future, it 
could induce the extraction of deeper and deeper flows, with a consequent 


increase in arsenic. 


The data obtained also allowed the determination of the Meteoric 
Water Line of the Cuernavaca Aquifer. This line will be very useful for 


future geohydrological and environmental studies. 


Acknowledgments 


The Mexican National Network of Isotopes in Precipitation (Known as 
RENIP) is a collaboration of the Mexican Institute of Water Technology, 
the National Meteorological Service, the Mexico City Water System, and 
the International Atomic Energy Agency, through project RLA7024 
"Integration of Isotopic Hydrology in National Water Resources 


Assessments”. 


References 


Brugge, D., € Buchner, V. (2011). Health effects of uranium: New 
research findings. Reviews on Environmental Health, 26(4), 231- 
24. 


Canter, L. W. (1997). Nitrate in groundwater. Boca Raton, USA: Lewis 
Publishers. 


Cardona, A. (2004). Salinization in coastal aquifers of arid zones: An 
example from Santo Domingo, Baja California Sur, Mexico. 


Environmental Geology, 45, 350-366. 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


pez Que: (Ml) Check for updates 
Tecnología y 
CienciaszAgua 
Clark, 1. D., 8 Fritz, P. (1997). Environmental isotopes in hydrogeology. 
New York, USA: CRE Press. DOI: 
https://doi.org/10.1201/9781482242911 


CDPH, California Department of Public Health. (2013). Drinking water 
contaminants: Nitrate. Sacramento, USA: California Department of 
Public Health. 


Conagua, Comisión Nacional del Agua. (2002). Determinación de la 
disponibilidad de agua en el acuífero Cuernavaca (1701), estado de 
Morelos. México, DF, México: Subdirección General Técnica, 


Gerencia de Aguas Subterráneas. 


Conagua, Comisión Nacional del Agua. (2010). Programa Hídrico Visión 
2030 del estado de Morelos. México, DF, México: Comisión Nacional 
del Agua. 


Conagua, Comisión Nacional del Agua. (2020). Actualización de la 
disponibilidad media anual de agua en el acuífero Cuernavaca 
(1701), estado de Morelos. México, DF, México: Subdirección 


General Técnica, Gerencia de Aguas Subterráneas. 


Dansgaard, W. (1964). Stable isotopes in precipitation. Te/llus, 16, 436- 
368. 


Del-Campo-Delgado, M. (2016). Lixiviación de agroquímicos en campos 
de golf bajo diferentes dotaciones de riego y sustratos (tesis 


doctoral). Universidad Autónoma del Estado de México, México. 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


pez Que: (Ml) Check for updates 
Tecnología y 
CienciaszAgua 
Esquivel-Hernández, G., Mosquera, G. M., Sánchez-Murillo, R., Quesada- 

Román, A., Birkel, C., Crespo, P., Célleri, R., Windhorst, D., Breuer, 
L., 8 Boll, J. (2019). Moisture transport and seasonal variations in 
the stable isotopic composition of rainfall in Central American and 
Andean Páramo during El Niño conditions (2015-2016). 
Hydrological Processes, 33(13), 1802-1817. 


Froehlich, K., Kralik, M., Rank, D., 8 Scheifinger, H. (2008). Deuterium 
excess in precipitation of Alpine regions. Moisture recycling. 


Isotopes in Environmental Health Studies, 44(1), 61-70. 


Galewsky, J., Steen-Larsen, H. C., Field R. D., Worden J., Risi C., 8 
Schneider M. (2016), Stable isotopes in atmospheric water vapor 
and applications to the hydrologic cycle. Reviews of Geophysics, 54, 
809-865. DOI: 10.1002/2015RG000512 


Gastmans, D., Santos, V., Aparecida-Galhardi, J., Felipe-Gromboni, J., 
Vianna-Batista L., Miotlinski, K., Kiang-Chang H., € Silvio-Govone, 
J, (2017). Controls over spatial and seasonal variations on isotopic 
composition of the precipitation along the central and eastern 
portion of Brazil. Isotopes in Environmental and Health Studies, 
53(5), 518-538. DOI: 10.1080/10256016.2017.1305376 


Gat, J. R. (2000). Atmospheric water balance. The isotopic perspective, 
Hydrological Processes, 14, 1357-1369. 


Gibbs, R. J. (1970). Mechanisms controlling world water chemistry. 
Science, 170(3962), 1088-1090. 


Open Access bajo la licencia CC BY-NC-SA 4.0 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


Tecnología y ciencias del agua, ISSN 2007-2422, 
15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 
CienciaszAgua 
Guo, Z., Yan, C., Wang, Z., Xu, F., € Yang, F. (2020). Quantitative 
identification of nitrate sources in a coastal peri-urban watershed 
using hydrogeochemical indicators and dual isotopes together with 


the statistical approaches. Chemosphere, 243, 125364. 


Haller, L., McCarthy, P., O'Brien, T., Riehle, J., 8 Stuhldreher, T. (2013). 
Nitrate pollution of groundwater. Oceanside, USA: Alpha Water 


Systems Inc. 


IAEA, International Atomic Energy Agency. (2002). A new device for 
montly rainfall sampling for GNIP. Water and Environment 
Newsletter, (16) special issue on the global network of isotopes in 


precipitation, 5. 


IMTA, Instituto Mexicano de la Tecnología del Agua. (2012). Plan integral 
para el manejo sustentable de las barrancas del norponiente del 
estado de Morelos. Jiutepec, México: Instituto Mexicano de la 


Tecnología del Agua, Fundación Gonzalo Río Arronte. 


INEGI, Instituto Nacional de Estadística y Geografía. (1991). Censo 
Agrícola, ganadero y ejidal del INEGI. Aguascalientes, México: 


Instituto Nacional de Estadística y Geografía. 


Jaimes-Palomera, L. R., Cortes-Silva, A., Vazquez-Sanchez, E., Aravena, 
R., Fritz, P., 8 Drimmie, R. (1989). Geoquímica isotópica del 
sistema hidrogeológico del valle de Cuernavaca, estado de Morelos, 
México. Geofísica Internacional, 28(2), 219-244. 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


open (oJaccess 0) Check for updates 
Tecnología y 
CienciaszAgua 
Kurttio, P., Komulainen, H., Leino, A., Salonen, L., Auvinen, A., 8 Saha, 
H. (2005). Bone as a possible target of chemical toxicity of natural 


uranium in drinking water. Environmental Health Perspectives, 113 
(1), 68-72. 


Liu, C. Q., Li, S. L., Lang, Y. C., 8 Xiao, H. Y. (2006). Using 9*N- and 
9180 values to identify nitrate sources in karst ground water, 
Guiyang, southwest China. Environmental Science € Technology, 
40(22), 6928-6933. 


Martinelli, L. A., Victoria, R. L., Sternberg, L. S., Ribeiro, A., € Moreira, 
M. Z. (1996). Using stable isotopes to determine sources of 
evaporated water to the atmosphere in the Amazon basin. Journal 
of Hydrology, 183, 191-204. 


Morales-Casique, E., Guinzberg-Belmont, J., € Ortega-Guerrero, A. 
(2016). Regional groundwater flow and geochemical evolution in 
the Amacuzac River Basin, Mexico. Hydrogeology Journal, 24, 1873- 
1890. 


Merlivat, L., € Jouzel, J. (1979), Global climate interpretation of the 
deuterium-oxygen 18 relationship for precipitation. Journal de 
Geophysical Research, 84(C8), 5029-5033. 


Nolan, J., € Weber, K. A. (2015). Natural uranium contamination in major 
U.S. aquifers linked to nitrate. Environmental Science and 
Technology Letters, 2(8), 215-220. 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


erez Que: (Ml) Check for updates 
Tecnología y 
CienciaszAgua 
Ortega, L. V., Garcia, S. A., € Pelayo, B. R. (2003). Aspectos 
geohidrológicos de los acuíferos del estado de Morelos. En: Oswald- 
Spring, U. El recurso agua en el Alto Balsas (pp. 93-107). México, 
DF, México: Universidad Nacional Autónoma de México, Centro 


Regional de Investigaciones Multidisciplinarias. 


Otte, 1., Detsch, F., Gútlein, A., Scholl, M., Kiese, R., Appelhans, T., 8 
Nauss, T. (2017). Seasonality of stable isotope composition of 
atmospheric water input at the southern slopes of Mt. Kilimanjaro, 
Tanzania. Hydrological Processes, 31(22), 3932-3947. DOI: 
10.1002/hyp.11311 


Pacheco, J., 8 Cabrera, A. (1997). Groundwater contamination by nitrates 


in the Yucatan Peninsula, Mexico. Hydrogeology Journal, 5, 47-53. 


Pacheco, J., Marín, L., Cabrera, A., Steinich, B., 8 Escolero, O. (2001). 
Nitrate temporal and spatial patterns in 12 water-supply Wells, 
Yucatán México. Environmental Geology, 40(6), 708-715. 


Paradis, C. J., Jagadamma, S., Watson, D. B., McKay, L. D., Hazen, T. C., 
Park, M., € Istok, J. D. (2016). In situ mobility of uranium in the 
presence of nitrate following sulfate-reducing conditions. Journal of 
Contaminant Hydrology, 187, 55-64. 


Pastén-Zapata, E., Ledesma-Ruiz, R., Harter, T., Ramírez, A. l., A 
Mahiknecht, J. (2014). Assessment of sources and fate of nitrate in 
shallow groundwater of an agricultural area by using a multi-tracer 
approach. Science of the Total Environment, 470-471, 855-864. 
DOI: 10.1016/j.scitotenv.2013.10.043 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 


(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 
CienciaszAgua 
Panno, S. V., Hackley, K. C., Hwang, H. H., Greenberg, S., Krapac, 1. G., 
Landsberger, S., € O'Kelly, D. J. (2002). Source identification of 
sodium and chloride contamination in natural waters. The 12th 


Annual Conference of the Illinois Groundwater Consortium, USA. 


Peng, H., Mayer, B., Harris, S., 8 Roy, K. R. (2007). The influence of 
below-cloud secondary effects on the stable isotope composition of 
hydrogen and oxygen in precipitation at Calgary, Alberta, Canada. 
Tellus B: Chemical and Physical Meteorology, 59(4), 698-704. 


POZCI, Programa de Ordenación de Zona Conurbada Intermunicipal. 
(2009). Programa de Ordenación de Zona  Conurbada 
Intermunicipal en su modalidad de Centro de Población. 
Cuernavaca, Emiliano Zapata, Jiutepec, Temixco y Xochitepec. 
Cuernavaca, México: Programa de Ordenación de Zona Conurbada 


Intermunicipal. 


Raisz, E. (1964). Landforms of Mexico. Mapa con texto, escala 1:3 000 
000 (2nd ed.). Cambridge, USA: Office of Naval Research. 


Ren, C., Zhang, Q., Wang, H., € Wang, Y. (2021). Identification of sources 
and transformations of nitrate in the intense human activity region 
of North China using a multi-isotope and Bayesian model. 
International Journal of Environmental Research and Public Health, 
18, 8642. DOI: https://doi.org/10.3390/ijerph18168642 


Ritter, W. F., € Chirnside, A. E. (1984). Impact of land use on ground- 
water quality in Southern Delaware. Ground Water, 22(1), 38-47. 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 


(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 
CienciaszAgua 
Sagarpa, Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y 
Alimentación. (2003). Evaluación de la alianza para el Campo 2002. 
Informe de evaluación estatal. México, DF, México: Secretaría de 


Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación. 


Secretaría de Desarrollo Sustentable de Morelos. (2009). Programa de 
Ordenación de Zona Conurbada Intermunicipal en su Modalidad de 
Centro de Población de Cuernavaca, Emiliano Zapata, Jiutepec, 
Temixco y Xochitepec. Cuernavaca, México: Secretaría de 


Desarrollo Sustentable de Morelos. 


SSA, Secretaría de Salubridad y Asistencia. (1994). NOM-127-SSA1- 
1994. Salud ambiental, agua para uso y consumo humano. Límites 
permisibles de calidad y tratamientos a que debe someterse el agua 
para su potabilización. México, DF, México: Secretaría de Salubridad 


y Asistencia. 


Singh, G., Sengor, S., Bhalla, A., Kumar, S. S., Stewart, B., Spycher, N., 
Ginn, T. M,, Peyton, B. M., 8 Sani, R. K. (2014). Reoxidation of 
biogenic reduced uranium: A challenge toward bioremediation. 
Critical Reviews in Environmental Science and Technology, 44(4), 
391-415. 


Suarez-Barragán, M. D., 8 Ureno-Luna, J. (1988). Calidad físico-química 
y bacteriológica del agua subterránea del Valle de Cuernavaca, 
Morelos (tesis de licenciatura). Facultad de Ingeniería, Universidad 


Nacional Autónoma de México, México. 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 
CienciaszAgua 
Scholl, M. A., Gingerich, S. B., € Tribble, G. W. (2002). The influence of 
microclimates and fog on stable isotope signatures used in the 
interpretation of regional hydrology: East Maui, Hawaii. Journal of 
Hydrology, 264, 170-184. DOI: https://doi.org/10.1016/S50022- 
1694(02)00073-2 


Stewart, M. K. (1975). Stable isotope fractionation due to evaporation 
and isotope exchange of falling water drops: Application to 
atmospheric processes and evaporation of lakes. Journal of 
Geophysical Research, 80, 1133-1146. 


Torres-Martínez, J. A., Abraham, M., Jurgen, M., Daesslé, L. W,, 
Cervantes-Avilés, P. A., 8 Ledesma-Ruiz, R. (2021). Estimation of 
nitrate pollution sources and transformations in groundwater of an 
intensive  livestock-agricultural area (Comarca Lagunera), 
combining major ¡ons, stable isotopes and MixSIAR model. 
Environmental Pollution, 269. DOI: 
https://doi.org/10.1016/j.envpol.2020.115445 


Victoria, R., Martinelli, L., Mortatti, J., 8 Richey, J. (1991). Mechanisms of 
water recycling in the Amazon Basin: Isotopic insights. Ambio: A 
Journal of Environment and Society, 20(8), 384-387. 


Viers, J. H., Liptzin, D., Rosenstock, T. S., Jensen, V. B., 8: Hollander, A. 
D. (2012). Nitrogen sources and loading to groundwater. 
Sacramento, USA: California State Water Resources Control Board. 

WHO, World Health Organization. (2007). Nitrate and nitrite in drinking- 


water. Background document for development of WHO guidelines 


for drinking-water quality. Geneva, Switzerland: World Health 


Organization. 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 1) Check for updates 
OPEN ACCESS 


Tecnología y 
CienciaszAgua 
WHO, World Health Organization. (2008). Guidelines for drinking-water 
quality (3"4 ed.). Geneva, Switzerland: World Health Organization. 


Zhao, Z. (2015). A global assessment of nitrate contamination in 
groundwater  (Internship report). Delft, The Netherlands: 


International Groundwater Resources Assessment Centre. 


2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 15(3), 186-249. DOI: 10.24850/j-tyca-2024-03-05 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


